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ABSTRACT. We consider a stochastic volatility stock price model in which the volatility is 
a non-centered continuous Gaussian process with arbitrary prescribed mean and covari¬ 
ance. By exhibiting a Karhunen-Loeve expansion for the integrated variance, and using 
sharp estimates of the density of a general second-chaos variable, we derive asymptotics 
for the stock price density and implied volatility in these models in the limit of large or 
small strikes. Our main result provides explicit expressions for the first three terms in 
the expansion of the implied volatility, based on three basic spectral-type statistics of the 
Gaussian process: the top eigenvalue of its covariance operator, the multiplicity of this 
eigenvalue, and the L 2 norm of the projection of the mean function on the top eigenspace. 
Strategies for using this expansion for calibration purposes are discussed. 
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1. Introduction 

In this article, we characterize the extreme-strike behavior of implied volatility curves 
for fixed maturity for uncorrelated Gaussian stochastic volatility models. This introduc¬ 
tion contains a careful description of the problem's background and of our motivations. 
Before going into details, we summarize some of the article's specificities; all terminology 
in the next two paragraphs is referenced, defined, and/or illustrated in the remainder of 
this introduction. 

We hold calibration of volatility smiles as a principal motivator. Cognizant of the fact 
that non-centered Gaussian volatility models can be designed in a flexible and parsimo¬ 
nious fashion, we adopt that class of models, imposing no further conditions on the mar¬ 
ginal distribution of the volatility process itself, beyond pathwise continuity. The spectral 
structure of second-Wiener chaos variables allows us to work at that level of generality. 
We find that the first three terms in the extreme-strike implied volatility asymptotics - 
which is typically amply sufficient in applications - can be determined explicitly thanks 
to three parameters characterizing the top of the spectral decomposition of the integrated 
variance. In order to prove such a precise statement while relying on a moderate amount 
of technicalities, we make use of the the simplifying assumption that the stochastic volatil¬ 
ity is assumed independent of the stock price's driving noise. 
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When considering the trade-off between this restriction and calibration considerations, 
we observe that our model flexibility combined with known explicit spectral expansions 
and numerical tools may allow practicioners to compute the said spectral parameters 
in a straightforward fashion based on smile features, while also allowing them to select 
their favorite Gaussian volatility model class. Specific examples of Gaussian volatility 
processes are non-centered Brownian motion. Brownian bridge, and Ornstein-Uhlenbeck 
models. This last sub-class can be particularly appealing to practitioners since it contains 
stationary volatilities, and includes the well-known Stein-Stein model. We also mention 
how any Gaussian model specification, including long-memory ones, can be handled, 
thanks to the numerical ability to determine its spectral elements. We understand that 
the assumption of the stochastic volatility model being uncorrelated implies the symme¬ 
try of the implied volatility in the wings, which in some applications, is not a desirable 
feature; on the other hand, in many option markets, liquidity considerations limit the abil¬ 
ity to calibrate using the large-strike wing (see the calibration study on SPX options in [20, 
Section 5.4]). The case of a correlated Gaussian stochastic volatility model is more com¬ 
plicated, but constitutes an interesting mathematical challenge, which we will investigate 
separately from this article, since one may need to develop completely new methods and 
techniques. An important step toward a better understanding of the asymptotic behavior 
of the implied volatility in some correlated stochastic volatility models is found in the 
articles [13, 14]. Another problem which is mathematically interesting and important in 
practice is the asymptotics for implied volatility in small or large time to maturity, on 
which we report in separate works. 

1.1. Background and heuristics. Studies in quantitative finance based on the Black-Scholes- 
Merton framework have shown awareness of the inadequacy of the constant volatility as¬ 
sumption, particularly after the crash of 1987, when practitioners began considering that 
extreme events were more likely than what a log-normal model will predict. Propositions 
to exploit this weakness in log-normal modeling systematically and quantitatively have 
grown ubiquitous to the point that implied volatility (IV), or the volatility level that mar¬ 
ket call option prices would imply if the Black-Scholes model were underlying, is now 
a bona fide and vigorous topic of investigation, both at the theoretical and practical level. 
The initial evidence against constant volatility simply came from observing that IV as a 
function of strike prices for liquid call options exhibited non-constance, typically illus¬ 
trated as a convex curve, often with a minimum near the money as for index options, 
hence the term Volatility smile'. 

Stock price models where the volatility is a stochastic process are known as stochastic 
volatility models; the term 'uncorrelated' is added to refer to the submodel class in which 
the volatility process is independent of the noise term driving the stock price. In a sense, 
the existence of the smile for any uncorrelated stochastic volatility model was first proved 
mathematically by Renault and Touzi in [30]. They established that the IV as a function of 
the strike price decreases on the interval where the call is in the money, increases on the 
interval where the call is out of the money, and attains its minimum where the call is at the 
money. Note that Renault and Touzi did not prove that the IV is locally convex near the 
money, but their work still established stochastic volatility models as a main model class 
for studying IV; these models continued steadily to provide inspiration for IV studies. 
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A current emphasis, which has become fertile mathematical ground, is on IV asymp¬ 
totics, such as large/small-strike, large-maturity, or small-time-to-maturity behaviors. 
These are helpful to understand and select models based on smile shapes. Several tech¬ 
niques are used to derive IV asymptotics. For instance, by exploiting a method of mo¬ 
ments and the representation of power payoffs as mixtures of a continuum of calls with 
varying strikes, in a rather model-free context, R. Lee proved in [28] that, for models with 
positive moment explosions, the squared IV's large strike behavior is of order the log- 
moneyness log (K/So) times a constant which depends explicitly on supremum of the 
order of finite moments. A similar result holds for models with negative moment explo¬ 
sions, where the squared IV behaves like K i-A log (So /K) for small values of K. More 
general formulas describing the asymptotic behavior of the IV in the 'wings' (K -A 0 or 
+oo) were obtained in [3, 4, 5,22, 23, 25, 17] (see also the book [ 11]). 

From the standpoint of modeling, one of the advantages of Lee's original result is the 
dependence of IV asymptotics merely on some simple statistics, namely as we mentioned, 
in the notation in [28], the maximal order p of finite moments for the underlying Sj, i.e. 


P(T) ■= sup 


p G 1R : E (Sj) 


\ P+1 


< oo 


This allows the author to draw appropriately strong conclusions about model calibration. 
A typical class in which p is positive and finite is that of Gaussian volatility models, which 
we introduce next. 

We consider the stock price model of the following form: 

dS t = rS t dt + \Xt\S t dW t : f e [ 0, T\, (1) 

where the short rate r is constant, X (f) = m(t) + X (t) with m an arbitrary continu¬ 
ous deterministic function on [0, T] (the mean function), and X is a continuous centered 
Gaussian process on [0, T] independent of W, with arbitrary covariance Q. Note that it is 
not assumed in (36) that the process X is a solution to a stochasic differential equation as 
is often assumed in classical stochastic volatility models. A well-known special example 
of a Gaussian volatility model is the Stein-Stein model introduced in [ 3], in which the 
volatility process X is the so-called mean-reverting Ornstein-Uhlenbeck process satisfy¬ 
ing 

dXt = a(m — Xf) dt + fidZt (2) 

where m is the level of mean reversion, oc is the mean-reversion rate, and ft is level of 
uncertainty on the volatility; here Z is another Brownian motion, which maybe correlated 
with W. In the present paper, we adopt an analytic technique, encountered for instance 
in the analysis of the uncorrelated Stein-Stein model by this paper's first author and E.M. 
Stein in [24] (see also [21]). 

Returning to the question of the value of p, for a Gaussian volatility model, it can some¬ 
times be determined by simple calculations, which we illustrate here with an elementary 
example. Assume S is a geometric Brownian motion with random volatility, i.e. a model 
as in (1) where (abusing notation) |Xf | is taken the non-time-dependent cr |X| where J is 
a constant and X is an independent unit-variance normal variate (not dependent on t). 
Thus, at time T, with zero discount rate, Sj = So exp (cr |X| Wj — a 2 X 2 T/2) . To simplify 
this example to the maximum, also assume that X is centered; using the independence of 
X and W, we get that we may replace | X| by X in this example, since this does not change 
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the law of Sj (i.e. in the uncorrelated case, X's non-positivity does not violate standard 
practice for volatility modeling). Then, using maturity T = 1, for any p > 0, the pth 
moment, via a simple change of variable, equals 

E [{S ' Y] = 2 ^TTW 2 //* ix iy exp H ( y2 + * ” 2 7fe"’ y )) 

which by an elementary computation is finite, and equal to (l + pa 1 ) / (l + pa 1 — p 2 a 2 ), 
if and only if 

~ , 1 /i r 

p<p + 1 -2 + Vi + ^2- 

In the cases where the random volatility model X above is non-centered and is correlated 
with W, a similar calculation can be performed, at the essentially trivial expenses of in¬ 
voking affine changes of variables, and the linear regression of one normal variate against 
another. 

The above example illustrates heuristically that, by Lee's moment formula, the com¬ 
putation of p might be the quickest path to obtain the leading term in the large-strike 
expansion of the IV, for more complex Gaussian volatility models, namely ones where 
the volatility X is time-dependent. However, computing p is not necessarily an easy task, 
and appears, perhaps surprisingly, to have been performed rarely. For the Stein-Stein 
model, the value of p can be computed using the sharp asymptotic formulas for the stock 
price density near zero and infinity, established in [24] for the uncorrelated Stein-Stein 
model, and in [14] for the correlated one. These two papers also provide asymptotic 
formulas with error estimates for the IV at extreme strikes in the Stein-Stein model. Be¬ 
yond the Stein-Stein model, little was known about the extreme strike asymptotics of 
general Gaussian stochastic volatility models. In the present paper, we extend the above- 
mentioned results from [24] and [14] to such models. 

1.2. Motivation and summary of main result. Adopting the perspective that an asymp¬ 
totic expansion for the IV can be helpful for model selection and calibration, our objective 
is to provide an expansion for the IV in a Gaussian volatility model relying on a minimal 
number of parameters, which can then be chosen to adjust to observed smiles. The re¬ 
striction of non-correlated volatility means that the stock price distribution is a mixture 
model of geometric Brownian motions with time-dependent volatilities, whose mixing 
density at time T is that of the square root of a variable in the second-chaos of a Wiener 
process. That second-chaos variable is none other than the integrated variance 

r T := f T X 2 ds. 

Jo 

By relying on a general Hilbert-space structure theorem which applies to the second 
Wiener chaos, we prove that, in the most general case of a non-centered Gaussian sto¬ 
chastic volatility with a possible degeneracy in the eigenstructure of the covariance Q of 
X viewed as a linear operator on L 2 ([0, T]) (i.e. when the top eigenvalue Aj is allowed 
to have a multiplicity n\ larger than 1), the large-strike IV asymptotics can be expressed 
with three terms which depend explicitly on T and on the following three parameters: Ai, 
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Hi, and the ratio 


S= \\P El m\\ /M 


where 11 Pe 1 m 11 is the norm in L 2 ([0, T ]) of the orthogonal projection of the mean function 
m on the first eigenspace of Q. Specifically, with I (K) the IV as a function of strike K, 
letting k := log (K/sq) — rT be the discounted log-moneyness, as k —* +oo, we prove 


I(K) = L(T,Ai)Vk + M (T, X\,5) - ^ + 0 (^=) (3) 

where the constants L and M depend explicitly on T and Ai, and M also depends explic¬ 
itly on 5. The details of these constants are in Theorem 4 on page 15. A similar asymptotic 
formula is obtained in the case where k -A — oo, using symmetry properties of uncorre¬ 
lated stochastic volatility models (see formula (50) on page 19). The specific case of the 
Stein-Stein model is expanded upon in some detail. Numerical illustrations of how the 
small strike asymptotics can be used are provided in Section 7, in the context of calibra¬ 
tion; we explain some of the calibration ideas in the next subsection. 


1.3. Practical implications. The first-order constant L is always strictly positive. The 
second-order term (the constant M) vanishes if and only if m is orthogonal to the first 
eigenspace of Q, which occurs for instance when m = 0. The third-order term vanishes 
if and only if the top eigenvalue has multiplicity one, which is typical. The behavior of L 
and M as functions of T is determined partly by how the top eigenvalue Ai depends on 
T, which can be non-trivial. In the present paper, we assume T is fixed. 

For fixed maturity T, assuming that Q has lead multiplicity n.\ = 1 for instance, a 
practitioner will have the possibility of determining a value Ai and a value 5 to match 
the specific root-log-moneyness behavior of small- or large-strike IV; moreover in that 


Jo ei(t)dt 


where 


y liquid at extreme 


case, choosing a constant mean function m, one obtains 5 = m 2 A-, 1 

e\ is the top eigenfunction of Q. Market prices may not be sufficient 
strikes to distinguish between more than two parameters; this is typical of calibration 
techniques for implied volatility curves for fixed maturity, such as the 'stochastic volatility 
inspired' (SVI) parametrization disseminated by J. Gatheral: see [18,19] (see also [20] and 
the references therein). Our result shows that Gaussian volatility models with non-zero 
mean are sufficient for this flexibility, and provide equivalent asymptotics irrespective of 
the precise mean function and covariance eigenstructure, since modulo the disappearance 
of the third-order term in the unit top multiplicity case n\ = 1, only A-| and 5 are relevant. 
Moreover, our Gaussian parametrization is free of arbitrage, since it is based on a semi¬ 
martingale model (1). In the case of SVI parametrization, the absence of arbitrage can be 
non-trivial, as discussed in [18] and [20]. 

Modelers wishing to stick to well-known classes of processes for X may then adjust 
the value of Ai by exploiting any available invariance properties for the desired class. For 
example, if X is standard Brownian motion, or the Brownian bridge, on [0, T], we have 
Ai = 4T 2 / n or Ai = T 2 / n respectively, and these values scale quadratically with respect 
to a multiplicative scaling constant for X, beyond which an arbitrary mean value m may 
be chosen. If X is the mean-zero stationary OU process, we have Ai = /3 2 / [coj + cc 2 ) 
where cvj is the smallest positive solution of 2amcos (cvT) + (a 2 — to) sin (coT) = 0, in 
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which case, for a fixed arbitrarily selected rate of mean reversion a., a scaling of Ai is 
then equivalent to selecting the variance of X, while a constant mean value m can then 
be selected independently. [9, Chapter 1] can be consulted for the eigenstructure of the 
covariance of Brownian motion and the Brownian bridge, which are classical results, and 
for a proof of the eigenstructure of the OU covariance (see also [ 1]). The top eigenfunc¬ 
tions in all three of these cases are known explicit trigonometric functions (see [9, Chapter 
1]), and need to be referenced when selecting m. For the OU bridge, the eigenstructure of 
Q (equivalently known as the Karhunen-Loeve expansion of Q) was found in [10], while 
in [12], such an expansion was characterized for special Gaussian processes generated by 
independent pairs of exponential random variables. On the other hand, fractional Brow¬ 
nian motion and OU processes driven by fractional Brownian motion do not fall in the 
class of Gaussian processes for which the Karhunen-Loeve expansion is known explic¬ 
itly. However, efficient numerical techniques exist to compute the eigenfunctions and 
eigenvalues in these cases: see [9, Chapter 2], OU processes driven by fractional Brow¬ 
nian motion were proposed early on for option pricing, and recently analyzed in [8, 7]. 
Section 7 illustrates how, in the case of the Stein-Stein model and of its version driven by 
fraction Brownian motion, the explicit, semi-explicit, or numerically accessible Karhunen- 
Loeve expansion of X can be used in conjunction with the asymptotics (3) for calibrating 
parameters. We find that market liquidity considerations limit the theoretical range of 
applicability of calibration strategies, but that significant practical results are nonetheless 
available. 

The remainder of this article is structured as follows. Section 2 sets up a convenient 
second-chaos representation for the model's integrated volatility. Section 3 uses calcula¬ 
tions in a proof in [6] to derive precise asymptotics for the density of the mixing distri¬ 
bution. Section 4 converts these asymptotics into two-sided estimates of the density of 
the stock price Sj, thanks to the analytic tools developed in [24, 21]. Sharper asymptotic 
formulas for the density can also be obtained, but they are not needed to derive sharp 
asymptotic formulas for the IV with three terms and an erorr estimate. In Section 5, we 
characterize the wing behavior of the implied volatility in Gaussian stochastic volatility 
models. The special case of the uncorrelated Stein-Stein model is studied in more detail in 
Section 6. Finally, our practical study of calibration strategies, with numerics, is in Section 
7. 

2. General setup and second-chaos expansion of the integrated variance 

Let X be an almost-surely continuous Gaussian process on a filtered complete probabil¬ 
ity space (O, J-, {X/ [, IP) with mean and covariance functions denoted by m(t) = E[X ; ] 
and 

Q(t,s) = cov{X t ,X s ) = E [(X, - m(f)) (X s - m(s))], 

respectively. While such processes used in a jump-free quantitative finance context for 
volatility modeling will require, in addition, that X be adapted to filtration of the Wiener 
process W driving the stock-price (as in (1)), under our simplifying assumption that X 
and W be independent, this adaptability assumption can be considered as automatically 
satisfied, or equivalently, as unnecessary, since the filtration of W can be augmented by 
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the natural filtration of X. Define the centered version of X : 

Xf := Xf — m{t), t > 0. 

Fix a time horizon T > 0. It is not hard to see that Q(s, s) > 0 for all s > 0. Since 
the Gaussian process X is almost surely continuous, the mean function t ha m(t ) is a 
continuous function on [0, T], and the covariance function (f, s) ha Q(t,s ) is a continuous 
function of two variables on [0, T] 1 2 . This is a consequence of the Dudley-Fernique theory 
of regularity, which also implies that m and Q boast moduli of continuity bounded above 
by the scale h ha log~ 1/2 ( h _1 ) (see [1]), but this can also be established by elementary 

l 

means . 

In our analysis, it will be convenient to refer to the Karhunen-Loeve expansion of 
X. Applying the classical Karhunen-Loeve theorem to X (see, e.g., [34], Section 26.1), 
we obtain the existence of a non-increasing sequence of non-negative summable reals 
{A w : n = 1,2,. .an i.i.d. sequence of standard normal variates {Z n : n = 1,2, ..and 
a sequence of functions {e n : n = 1,2,...} which form an orthonormal system in L 2 ([0, T]), 
such that 

OO _ 

Xf = X] V%ie n (t)Z n . (4) 

n—1 

In (4), {e n = e n j} are the eigenfunctions of the covariance Q acting on L 2 ([0, T]) as the 
following operator 

T 

£(/)(*)=/ f(s)Q(t,s)ds, /GL 2 ([ 0,r]), 0 < f < T, 

Jo 

and { A n = A n j }, n > 1, are the corresponding eigenvalues (counting the multiplicities). 
We always assume that the orthonormal system {e n } is rearranged so that 

M = A 2 — • • • = A Wl > A Wl _)_2 = ^ni+2 = • • • = A„ 1+ „ 2 > ... 

In particular, Ai is the top eigenvalue, and n\ is its multiplicity. 

Using (4), we obtain 

rT _ pT ( 00 _ \ oo 

/ Xfdt= E yfKe n {t)Zn ) dt = J] AnZ 2 . (5) 

■ y 0 V«=l / n=! 

It is worth pointing out that this expression for the integrated variance of the centered 
volatility Jq Xjdt is in fact the most general form of a random variable in the second 
Wiener chaos, with mean adjusted to ensure almost-sure positivity of the integrated vari¬ 
ance. This is established using a classical structure theorem on separable Flilbert spaces, 
as explained in [29, Section 2.7.4]. In other words (also see [29, Section 2.7.3] for additional 

1 

The continuity of the process X implies its continuity in probability on O. Hence, the process X is 
continuous in the mean-square sense (see, e.g., [27], Lemma 1 on p. 5, or invoke the equivalence of LT 
norms on Wiener chaos, see [29]). Mean-square continuity of X implies the continuity of the mean function 
on [0, T], In addition, the autocorrelation function of the process X, that is, the function R(t,s ) = E [XfX s ], 
(t,s) G [0, T] 2 , is continuous (see, e.g., [2], Lemma 4.2). Finally, since Q(t,s) — R(t,s ) — m(t)m(s), the 
covariance function Q is continuous on [0, T] 2 . 
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details), any prescribed mean-adjusted integrated variance in the second chaos is of the 
form 

^ (T) ■■= JJ G (M) dZ (s) dZ (f) + 2 II G|| 2 t2(laTp) 

for some standard Wiener process Z and some function G £ L 2 ([0, T] 2 ), and moreover 

one can find a centered Gaussian process X such that V (T) = Xjdt and one can com¬ 
pute the coefficients A„ in the Karhunen-Loeve representation 5 as the eigenvalues of 
the covariance of X. When using the non-centered process X, this analysis immediately 
yields 



Xfdt 


J ^Xt + m(t )^ dt 
00 00 

E A » z » + 2 E 

n=l n—1 

+ r - 


m(t)e n (t)dt 


Z r 


/ o 


m(t) l dt. 


Set 

5 n = d n T = / m(t)e n (t)dt, n> 1, 
Jo 

and 

T 

s = s(T) = / m(t) 2 dt. 

Jo 

Then Bessel's inequality implies that 


( 6 ) 

(7) 

( 8 ) 


00 


E S n<S- 


(9) 


Denote 

00 

T = s - E 3 n- ( 10 ) 

n—1 

It is not hard to see that if the function t i—)• m(t ) belongs to the image space /C(L 2 [0, T]), 
then 

OO 

E s l = (ii) 

n—1 

and hence, T = 0. For instance, the previous equality holds for a centered Gaussian 
process X. 

Equality (6) can be rewritten as follows: 





n=1 



+ S 


E A > 


n=1 


z n + 




( 12 ) 
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It follows from (11) and (12) that if the function t > m(t) belongs to the image space 
JC(L 2 [0, T]), then 


T 


Xjdt 


00 


E A » 

n —1 



(13) 


Let us denote the noncentral chi-square distribution with the number of degrees of 
freedom k and the parameter of noncentrality A by y 2 (/c, A) (more information on such 
distributions can be found in [ 1] or in any probability textbook). Define a random vari¬ 
able Z n by 


Z 


n 


Z n + 



~ e2 

It is clear that Z n is distributed as ^ 2 (1, ^). Set 

P 1 = Zl, P2 = A ;1l _|_2/ P3 = +172 + 1/ • • • / 

and 


A r = — ( / X 2 dt - r 

Ai \J o 


Then, using (12), we see that 


At =x2{ni 'iP") + *z 2 ( n ^Xj ]+ ' 


(14) 


(15) 


where the repeated chi-squared notation is used abusively to denote independent chi- 
squared random variables. 


3. Asymptotics of the mixing density 

The asymptotic behavior of the complementary distribution function of an infinite lin¬ 
ear combination of independent central chi-square random variables was characterized 
by Zolotarev in [35]. Zolotarev's results were generalized to the case of noncentral chi- 
square variables by Beran (see [6]). Beran used some ideas from [26] in his work. We will 
employ Beran's results in the present paper. More precisely. Theorem 2 in [ 6 ] will be used 
below. This theorem provides an asymptotic formula for the complementary distribution 
function of an infinite linear combination of independent noncentral chi-square random 
variables. A sharper formula for the distribution density cjj of such a random variable can 
be extracted from the proof of Theorem 2 in [6] (see the very end of that proof). Adapting 
Beran's result to our case and taking into account estimate (9) (this estimate is needed to 
check the validity of the conditions in Beran's theorem), we see that 





A 


(16) 
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as x —E oo. In (16), the number A is given by 

Ai 


^ = n 

k=2 

x exp 


M ~ pk 

-1 OO 

E 


2 


1 


2 fc = 2 Ai Pk \ n= 


«iH-E n k 

E 


n lH- 


Set 


Then (16) gives 


«1 


' = ^4 


1 n=l 


<7r(*) = Ap^ (x;ni,<S) (l + 0 (x . 

For A > 0, the following formula is known: 

P x 2 ( x ' n '^) = \{j) i 2 *"*^*5-1 ("'/A*)/ * > 0 , 

where Z y is the modified Bessel function of the first kind (see, e.g., [21], Theorem 1.31). It 
is easy to see that formula ( 20 ) and the formula 


(17) 

(18) 

(19) 

( 20 ) 


7i/(f) — 


\/1nt 


l+O t 


i-i 


t —* oo. 


describing the asymptotic behavior of the I-Bessel function, imply that 


p x2 (x;n, A) = 


1 


2y/2n 

as x —> oo. On the other hand, if A = 0, then 


ft-1 )i-3 ,/i r _x±A . _ , _* 

A 4 x 4 e VA e 2 (l + O x 2 


p x2 (x;n,0) = 


1 

—x 2 exp 


x > 0 , 


( 21 ) 


( 22 ) 


(see, e.g.. Lemma 1.27 in [21]). 

Recall that we denoted by qj the distribution density of the random variable A 7 defined 
by (14). It follows from (19) and (21) that 


, \ A r _Nzl !+ JTZ _*±<E 
qr{x) = —^<5 4 x 4 e v e 2 


l + O X “2 


(23) 


2V2n 

as x -A 00 . The constants A and S in (23) are defined by (17) and (18), respectively For a 
centered Gaussian process X, (19) and (22) imply that 


q T (x) = 


A 


m-2 


2 ?r(a) 


X 


x 2 exp -- l + O x 2 


x > 0 , 


as x -A 00 , where 


^ = n(+^- 

k=2 \A-l~Pk 


(24) 


( 25 ) 


Indeed, in this case, we have s = 0, 8 n = 0 for all n > 1, £ = 0, and r = 0. 
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Our main goal is to characterize the asymptotic behavior of the distribution density pj 
of the random variable 

r T 

T t = / xfdt. ( 26 ) 

Jo 

It follows from (10) and (14) that T^ = + r. Therefore, 

Pt(x ) = (x^ x - T )) • ( 27 ) 

Theorem 1. Let pj be the distribution density of the random variable Y j given by (26). If 5 > 0, 
then the following asymptotic formula holds: 


p T (x) = Cx"'< 3 exp | exp { -^ 


x (1 + O (x 2 


as x — > oo, where 


C = 


“I 


A 


2y/2n 


V 


1 / ”i 

E*5 

n —1 


exp 


v^oo r2 

s l—m=1 1 


2Ai 


(28) 


(29) 


In (28) and (29), f/ze constants 5 n , s, A, and S are defined by (7), (8), (17), and (18), respectively. 
On the other hand, for a centered Gaussian process X, we have 


Hi- 2 


as x —* 00 , ze/zere 


Pt(x) = Cx^““ exp 


C = 


x 


2 Ai 


1 + 0 X “2 


A -N 

A 2 
A l 


2 ^ r (+) 


Proof. Formulas (23) and (27) imply that 


Pt(x) = 


A 1 a 




n l~ l 


2\/2n Ai 


A x 4 <T^ 


Hi -1 


X 


n 1 

E^ 

v«=l 


Hi —3 


A -* 


1 ex p 


r — +” 1 A 2 


ni —3 


x (x — r) 4 exp 


+(x — r) 


2Ai 


exp 


x ) 

2 Ai / 


(30) 


(31) 


x 1 + O x 2 


as x -A 00 . 

Next, taking into account the formulas 


n 1 —3 —3 


(x —r) ±4 = x " 4 (l + 0(x x )) 


( 32 ) 
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and 

exp =exp { (1+0(I_i)) ' 
and simplifying the expression on the right-hand side of (32), we obtain formula (28). The 
proof of formula (30) is similar. Here we use (24) and (27). □ 

The next assertion follows form Theorem 1. 


Corollary 2. Let pj be the distribution density of the random variable T j given by (26). Then the 
following are true: 

(1) Ifn\ = l, then 


Pt(x) =Cx 


1 

2 


exp 



x 


( 1+0 (* 4 )) 


x 

2Ai 


as x — > oo. 

(2) Suppose X is a centered Gaussian process with ry = 1. Then 
p T (x) = Cx~hx p (l + 0 

as x -A oo. 


(33) 


(34) 


4. Stock price asymptotics 

In this section, we study stochastic volatility models, for which the volatility is de¬ 
scribed by the absolute value of a Gaussian process. 

Recall that in the present paper we assume that the asset price process S satisfies the 
following linear stochastic differential equation: 

dSt = rStdt + \Xt\SfdWt, (35) 

where X is a continuous Gaussian process on (O, JF, {J-)},P). In (35), W is a standard 
Brownian motion on (O, LF, P) with respect to the filtration { Tt }, and the symbol r stands 
for the constant interest rate. It will be assumed that the processes X and W are indepen¬ 
dent. The initial price of the asset will be denoted by so- 
Since (35) is a linear stochastic differential equation, we have 

St = so exp jrf — ^ J Xgds + J |X s |dW s 

The previous equality follows from the Doleans-Dade formula (see [ ]). Therefore, the 

discounted asset price process is given by the following stochastic exponential: 

S t = e~ rt S t = s 0 exp j ~ X 2 s ds + J* |X s |dW s j . (36) 

The next assertion states that under the restrictions on the volatility process used in the 
present paper, we are in a risk-neutral environment. 
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Lemma 3. Under the restrictions on the volatility process X in (35), the discounted asset price 
process S is a {(Ft}-martingale. 

Proof. Fix a time horizon T > 0. It suffices to prove that there exists 5 > 0 such that 

L = sup E[exp{<5X^}] < oo (37) 

o <t<T 

(see, e.g., [21], Corollary 2.11). It will be shown that (37) holds provided that 

1 


5 < 


2 max 0 <f<T Q(t, t) 


(38) 


We have 


L = _ sup 


< 


V2n 0<t<T \/Q{t,t) j ]R 

. cxn f - (()2 

1 eX P \ 2 Q(t,t) 

._ sup - , 

V2no<t<T y/Q(t,t) 


exp by 1 - 


R 

_ m(t) 

QM 

m(t) 2 


exp 


2Q(t,t) 
1 


2 Q(t,t) 


(y-m(t)f 


~ 5 )y 2 + PFruV \ d v- 


Q(t,ty 


Set a(t) = ~ 5 and /3(f) = ffjf- Then, completing the square in the previous 

exponential, we obtain 


2 Q(f,f) 


1 ex P 

L < __ sup - ,_ 

V^n o<t<r 


exp 


f Pi*) 2 ) 
IMOJ 


/ R 


exp C«W| 


V5T oSIt vWJ) l Mt) 2 Q(t, t) 

(3(t) 2 m(t ) 2 


exp 


dy 

(f(t) 2 m{t) 2 


I R 


exp | -a(f)y 2 | dy 


= a\ sup 


= #2 


o <t<T \iQ(t,t)oc(t ) 

1 


exp 


4a(f) 2Q(f, t) 

f m(t) 2 r 


0<KT \/l — 2SQ(t, t) 6XP 12Q(f, f) |_1 — 23Q(t, t) 


- 1 


= SUp 

< fl2 


1 


o<t<T \/l - 2 SQ(t,t) 

1 


exp 


bm(ty 


1 — 28Q(t, t) 
exp 


bm(ty 


v/W^I^maxo^rQlhT) 11 — 2^maxo<t<T 

Finally, using (38) and the estimate |m(f) | < M for all 0 < t < T, we obtain (37). 
This completes the proof of Lemma 3. 


□ 
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Since the processes X and W are independent, the following formula holds for the 
distribution density D t of the asset price Sp 

D t (x) = 


\/2nt 


1°° 1 /—1 

[ 

lo y exp 1 



l o 8 2 g if 

2fy2 8 


Pt(y)dy. 


(39) 


In (39), pt is the distribution density of the random variable 


Yt = <- 


1 f* 


t Jo 


Xfds 


The function pt is called the mixing density The proof of formula (39) can be found in 
[ 21] (see (3.5) in [21]). 


It is not hard to see that 


pt(y) = 2t ypt (V) , 

(40) 

where the symbol pt stands for the density of the realized volatility Yf 
follows from formula (28) that 

= /o' X?rfs. It 

My) = AyM~ exp {By} exp { -Cy 2 } 


x (l + o(y->)) 

(41) 

as y —> oo, where 


~ «i+i ~ / St ~ t 

^ = 2Ct ^' B = Vv c = 2v 

(42) 


Our next goal is to estimate the function Df. The asymptotic behavior as x —> oo of the 
integral appearing in (39) was studied in [ 4] (see also Section 5.3 in [24]). It is explained 
in [24] how to get an asymptotic formula for the integral in (39) in the case where an as¬ 
ymptotic formula for the mixing density is similar to formula (41). We refer the reader 
to the derivation of Theorem 6.1 in [21], which is based on formula (5.133) in Section 5.6 
of [21] and Theorem 5.5 in [21]. The latter theorem concerns the asymptotic behavior of 
integrals with lognormal kernels. Having obtained an asymptotic formula for the distri¬ 
bution density of the asset price, we can find a similar asymptotic formula for the call 
pricing function C at large strikes, and then obtain an asymptotic formula for the implied 
volatility I (see Section 10.5 in [21]). However, the coefficients in the fourth term (and 
the higher order terms) in the asymptotic expansion of the implied volatility I are com¬ 
plicated and not suitable for computations. For that reason, we will restrict ourselves to 
asymptotic formulas for the implied volatility with three terms and an error estimate. To 
obtain such formulas, it is not necessary to use sharp asymptotic for Df or C. It suffices to 
establish two-sided estimates for Df and C (see [21], Chapters 9 and 10). The discussion 
below exploits the following two-sided estimate for the mixing density pt (this estimate 
follows from (41): 

My) ~ exp |By| exp | -Cy 2 } 


(43) 
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as y -A oo. Using formula (39) and reasoning as was described above, we see that for 
T > 0, 

D TW «('ogH_) 4 ex P { r j r) 1 } 

_(3 | Vs£+t \ 

x v 2 ^ ) (44) 


as x -A oo. In (44), 



(45) 


The proof of (44) uses (43) and is similar to the proof of formula (6.49) on pages 186-187 
of [21]. 


5. Asymptotics of the implied volatility 


The call pricing function in the stochastic volatility model described by (35) will be 
denoted by C. We have 

C(T,K) =e~ rt E[(S T -K) + ] 

where T is the maturity and K is the strike price. We will fix T, and consider C as the 
function K ha C(K) of only the strike price K. The Black-Scholes implied volatility asso¬ 
ciated with the pricing function C will be denoted by I. More information on the implied 
volatility can be found in [19, 21]. 

The next assertion provides an asymptotic formula for the implied volatility in the sto¬ 
chastic volatility model given by (35). 


Theorem 4. The following formula holds for the implied volatility K ha I(K) as K -A oo; 


m = f 

T 4 


+ 


1 

,_\ 2 


8C + T+ VtY - ( y 8C + T - \/T j 


1- 
\ 2 


' lQ g U§T 


s 0 e r 


Y2B 


T2(8C + T)i 


YsC + T - v/T ) 2 ( YsC + T + Yt 


+ ^lo^ +o(fiog ^ 


lo S^ 


4 K 

The constants in (46) are defined by 

B = 


s 0 e' 


T n l _ 

- ~ y 5}, and C 

(Ai ftl 


T 

2M’ 


(46) 


(47) 


Proof The equalities in (47) follow from (18) and (45). 

We will use the following assertion which was established in [ 1] (see Theorem 10.17 
in [ 21 ]). 
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Theorem 5. Suppose that the density of the asset price Dj satisfies 


Dt(x) ~ x K h(x) 


as x —» oo, where a. < —2 and h is a slowly varying function. Then 


\/2 / 11 1 
m = 7f f° SK + lOS K^^)-T OS ^K^h(K) 


v/2 W - ^iogiog 1 


y/f y 5 K*+ 2 h(K) 2 

+ o((logK)" 


K a+2 h(K) 


as K —> oo. 


Set a. = — 


\/8C+T 

2 VT 


2 T ~~—— ] and 


Mi —3 


H x ) = ( log ) exp <| A /log * 


s 0 e 


T4(8C + T); 


s 0 e 


rT 


It is easy to see that a. < — 2 and h is a slowly varying function. Next, taking into account 
(44) and applying Theorem 5, we obtain 


,{K) = 7? V ^ ~ ^ 108 ^+ 108 mV i 108108 

■ ^ V ' 2> IOg ^ + l0g ft(^) _ 1 IOgl ° S ^f4(K) 


+ °l ( lo S^T 
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as K —» oo. Now, using the mean value theorem, we get 


T , r . V2 I, K , 1 1 , K 

W = -W (-* -1) lQ g + lo § - 2 lo § lo s 


x/T 


Sqc 


s 0 p' 


■ ^ V ( “ a “ 2) log ^ +log “ ^ log log ^ 


+ 0 log 




s 0 e 


rT 


lo g ]/ 1 + S l( K ) 


x/2 
' x/T 

- 2^108 ^ \/l+ S2 (K) 

+ °(( l0g ^ 


as K —> oo, where 


Si(KC) 


lQ g^)-I lQ g lQ g^T 

(-«-l) lo g^ 


and 


s 2 (K) 


(-*-2) log Jjr 


(48) 
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Next, using the formula \/l + s = 1 + + 0(s 2 ) as s —» 0 in (48), we obtain 


y /2 r _■ 

I(K) = v 7 -a - 1 - v 7 -a - 2 


Vr l 


1 

+ -^= 


ll °g^T 


1 


1 


s 0 e' 

lo S 777XT - 2 lo § lo S 


\JTT _ \J—a — 1 — 

K 


lo Siih 


+ °l ( log s 0 ^ 


V2 

VT 


+ 


V ~0i — 1 — \ f —Oi — 2 

1 1 


1108 ^ 


+ 


V 7 — Oi — 2 y/—a — 1 
n l 


B 


ri(8C + r)i 


l-«, lo S log ^7 


l0 S^T 


+ O log 


K 


s 0 e 


l-T 


as K -A oo. Finally, taking into account the definition of the parameter a, we see that (46) 
holds. 

The proof of Theorem 4 is thus completed. □ 


In terms of the log-moneyness k = log 


Theorem 4 can be formulated as follows. 


Theorem 6. The following formula holds for the implied volatility k i—>■ I(k) as k —»■ oo; 


m = f 

T 4 


+ 


8C + T+ - ^\J 8C + T — Vf^j 2 


Vk 


V2B 


+ 


T2(8C + T)4 
1 — Mi log k 


Vsc + T - s/t) 1 ( Vsc + T + \/r X 2 


Vk 


k + o n 


(49) 


Note that if Mi = 1, then the third term in the asymptotic expansion of the implied 
volatility in formula (46) vanishes. On the other hand, if the Gaussian process X is cen¬ 
tered, then B = 0, and the second term on the tight-hand side of (46) vanishes. The next 
corollary takes those facts into account. 
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Corollary 7. The following are true: 
(i) then as K -A oo. 


m = -tt 


2Af 


T 2 

+ 


(\/4 + Ai + \/Ai) 2 + (-\/4 + Ax — \/4 t) 2 

y/25 


' lo s^f 


s 0 e' 


T 2 (4 + Ai) 4 (y / 4+ _ Ai + \/Ai) 2 + (\/4 + Ai — \/Ai) 2 
K 


+ ° I I lo 8 


(ii) If X is a centered Gaussian process, then as K -A oo, 

W = —pr 


2Af 


T 2 


(\/4 + Ai + 1 /A 1) 2 + (\/4 + Ai — \/Ai) 2 

K 


+ 


i-m lo giog^r 


lo § Jr 


+ 0 log 


X 


s 0 e 


,rT 


(iii) If X is a centered Gaussian process and n\ = 1, then as K -A oo, 

W = T7 


2A? 


(\/4 + Ai + \/Ai) 2 + (\/4 + Ai — \/Ai) 2 
K 


Ti 

+ 01 (log 


s 0 c 


rT 



K 

s 0 e rT 



K 

s 0 e rT 


The constants in Corollary 7 are obtained from those in Theorem 4, by taking into ac¬ 
count (47). These constants depend only on the following characteristics of the process 
X: The largest eigenvalue Ai of the covariance operator fC on [0, T], and the Fourier coef¬ 
ficient 5\ of the mean function m with respect to the corresponding eigenfunction e\. We 
refer to Section 1.3 for further discussion. 

Since the processes X and W in (35) are independent, the stochastic volatility model 
described in (35) belongs to the class of the so-called symmetric models (see Section 9.8 in 
[21]). It is known that for symmetric models. 


m 



for all K > 0. 


(50) 


The next statements can be established using Theorem 4, Corollary 7, and formula (50). 
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Theorem 8. The following formula holds for the function K ha I(K) as K -A 0: 


m 


+ 


\ 2 


8 C + T+ Vt\ - U/8C + T- Vt\ 


\ 2 


log 


s 0 e 


rT 


K 


V2B 


+ 


T3(8C + T)? 


1 — fti log log 


\/ 8 C + T - Vt) 1 ( VSC + T + VT 


i 

/=\ 2 


spe 1 


.rT 


K 


soe 


.rT 


+ 0 log 


s 0 e 


rT 


lo g^ 

The constants in (51) are the same as in Theorem 4. 

Corollary 9. The following are true: 

(i) lfn\ — l, then as K -A 0, 


K 


m = - 


2A? 


T 2 

+ 


(\/4 + Ai + \/Ai )2 4 - (-\/4 + Ai — \/Ai )2 

V2h 


!og 


s 0 e 


rT 


K 


T2(4+Ai)5 


+ 0 log 


s 0 e 


rT 




(ii) IfX is a centered Gaussian process, then as K ^ 0, 

W = -TT 


2A? 


T 2 


(\/4 + Ai + i/Ax ) 2 + (\/4 + Ai — \/Ai ) 2 


!og 


S 0 T 


rT 




+ 


1 - ni log log 


a rT 


K 


Spg' 


rT 


+ 0 log 


S 0 G 


fC 


lo gTT 

(iii) If X is a centered Gaussian process and n\ = 1, then as K -A 0, 


W = tt 


2 a) 


T2 


(\/4 + Ax + \/Ai ) 2 + (\/4 + Ax — \/Ax ) 2 


!og 


s 0 e 


.rT 


K 


+ 0 log 


(51) 


(52) 


soe 


rT 


K 


6. Implied volatility in the uncorrelated Stein-Stein model 

The classical Stein-Stein model is an important special example of a Gaussian stochas¬ 
tic volatility model. The Stein-Stein model was introduced in [ 33 ], The volatility in the 
uncorrelated Stein-Stein model is the absolute value of an Ornstein-Uhlenbeck process 
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with a constant initial condition itiq. In this section, we also consider a generalization of 
the Stein-Stein model, in which the initial condition for the volatility process is a random 
variable Xo. Of our interest in the present section is a Gaussian stochastic volatility model 
with the process X satisfying the equation dXt = q(m — Xt)dt + crdZt. Here q > 0, m > 0, 
and cr > 0. It will be assumed that the initial condition Xo is a Gaussian random variable 
with mean mo and variance Oq, independent of the process Z. It is known that 

X f = e~ qt Xo + (1 - e~ qt )m + ae~ qt t e qu dZ U/ t > 0. (53) 

Jo 


If ctq = 0, then the initial condition is equal to the constant rtiQ. The mean function of the 
process X is given by 

m(t) = e~ qt m.Q + (1 — e~ qt )m, (54) 

and its covariance function is as follows: 


Q(t,s) = e~ q ^ j( 7 0 2 + ^ (>mm(Ls) _ ^ J . 


Therefore, the following formula holds for the variance function: 

r 2 / J2 ' 


cr" 

2 q 


cr} = — + e~ 2qt 


Vo ~ ~ 


2 q 


and hence, if <7 q = then the process X f — m(t), t G [0, T], is centered and stationary In 
this case, the covariance function is given by 


Q (t,s) = 

The Karhunen-Loeve expansion of the Ornstein-Uhlenbeck process is known explicitly 
(see [11]). Denote by w n the increasingly sorted sequence of the positive solutions to the 
equation 

a 2 w cos(zpT) + (qcr 2 — w 2 Oq — q 2 cr { g) sin(wT) = 0. (55) 

If Go = 0, then the equation in (55) becomes 

w cos (wT) + q sin (wT) = 0. (56) 


For the OU process in (53) with cr 0 ^ 0, we have n k = 1 for all k > 1; 


Ah — 


<7 


u>l\ + q 2 

for all n > 1; and 

e n (t) = K n [oQW n cos(w„f) + (cr 2 — qa { 2 ) sin(w„f)] 
for all n > 1 and t G [0, T], The constant K n in (58) is determined from 


1 


in = yrr^oV 2 - qv'o)(l - cos(2w„r)) + > 0 4 w^ ( T + y£-sin(2 w„T) 


,4_„2 


K 2 2 w 


2w n 


+ y(^ ~JVq) ( T- — sin(2w„r) 


(57) 

(58) 


(59) 
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for all n > 1. On the other hand, if gq = 0, then A„ is given by (57), while the functions e n 
are defined by 

1 

e n (t) = — - sin (w n t) (60) 


sm(2w n T) 


for all n > 1 and t E [0, T]. 

By the Karhunen-Loeve theorem, the Ornstein-Uhlenbeck process X in (53) can be rep¬ 
resented as follows: 

00 _ 

X t = e~ qt m 0 + (1 - e~ qt )m + ^/X^e n (t)Z n 

n —1 

where {Z n } n >\ is an i.i.d. sequence of standard normal variables. The eigenvalues A fJ , 
n > 1, and the eigenfunctions e n , n > 1, are given by (57) and (58) if Go 7 ^ 0, and by (57) 
and (60) if Go = 0. Recall that the numbers w n , n > 1, in (57) are solutions to the equation 
in (55) if Go f 0, and to the equation in (56) if Go = 0. We refer the interested reader to [ 1] 
for more details. 

Our next goal is to discuss the constants in the asymptotic formulas for the implied 
volatility at extreme strikes in the Stein-Stein model (see the formulas in Corollaries 7 
and 9). Since n\ = 1 for any OU process, the only parameters needed to compute the 
above-mentioned constants are Ai and 5\. If Gq 7 ^ 0, then we have 


Ai = 


G 


w\ + q 1 ' 


(61) 


where W\ is the smallest strictly positive solution to the equation in (55). 

The next assertion provides explicit formulas for the number 5\ = | () m(t)e-\ (t)dt. 

Lemma 10. (i) For the generalized nncorrelated Stein-Stein model with gq 7 ^ 0, 

Kim(<j 2 — qo~Q ) (1 — cos(wiT)) 


h = 


W\ 


+ Kicr 2 (mQ — m ) 


Wi[l 


- p-qT 


+ XiGq sin(ifiT) [(mo — m)e ^ + m] 
cos(wiT)] — qe-^ 7 sin(wiT) 


q 2 + w\ 


(62) 


where the constant K\ is determined from (59) with n = 1. The symbol w\ in (62) stands for the 
smallest strictly positive solution to (55). 


(ii) For the uncorrelated Stein-Stein model with Xq = mo F -almost surely, 

mq 2 ( 1 — cos(wiT)) + w\(m.Q — m cos(wiT)) 


h = 




Proof Taking into account (54) and (58), we see that 


8\ — b\ 


f T 

cos(wit)dt + &2 / e qt cos(wit)dt 

Jo 

[T fT 

+ &3 / sin(wit)dt + &4 / e qt sm(wit)dt, 
Jo Jo 


(63) 


(64) 
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where 

b\ = mK\OQW\ r b 2 = (mo - m)KiOoWi, 

&3 = mKiio 2 — q&o), and b 4 = (rag — m)Ki(a 2 - qa^). 
It remains to evaluate the integrals in (64). We have 

sin(zPiT) 


/ 0 


cos(w\t)dt = 


W\ 


' T e-*cos (w lt )dt = ^-^cos^nl+^-^sin^T) 


'0 


and 


^ 2 + w 2 

f T sm(w 1 t)dt= 1 - C ° s( - W ' T \ 

1 0 Wi 


" T sm(w\t)dt = Wl[1 - ~ 

3 <7 2 + wf 


(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 


In the proof of (67) and (69), we use the integration by parts formula twice. Now, taking 
into account formulas (64-69) and making simplifications, we establish formula (62). 
Next, suppose Cq = 0. Then (62) implies that 


61 = 


m 


1 — cos(zpiT) 


+ 


sin(2z<’iT) 10\ 

Wi[l — e -t?T cos(wiT) — c\e~^ T sin(zPiT) 


4 u>i 

mo — m 


sin(2w 1 T) 
4 w\ 


q 2 + W\ 


where w\ denotes the smallest strictly positive solution to (56). It is not hard to see, using 
the equality W\ cos (w\T) + q sin(zPiT) =0, that 


h = 


m 


1 — cos(zpiT) 


sin(2wiT) 
4 w\ 


W\ 


+ 


mo — m 


W\ 


fr sin(2zyiT) q 2 + w 2 ’ 

V 2 4wi 


and it is clear that (70) and (63) are equivalent. 
This completes the proof of Lemma 10. 


(70) 


□ 


Remark 11. Since for the generalized Stein-Stein model with a random initial condition we have 
n\ = 1, one can use the asymptotic formulas in the first parts of Corollaries 7 and 9 to characterize 
the wing-behavior of the implied volatility. The dependence of the parameters Ai and 5\, appearing 
in those formulas, on the model parameters is described in (61), (62), and (63). Originally, sharp 
asymptotic formulas for the implied volatility at extreme strikes in the uncorrelated Stein-Stein 
model with Xo = mo were obtained in [ i-] (see also Section 10.5 in [ 1]). However, explicit 
expressions, obtained in [24, 21] for the coefficients in the asymptotic formulas for the implied 
volatility in the Stein-Stein model, are significantly more complicated than those found in the 
present paper. 


















24 


ARCHIL GULISASHVILI, FREDERI VIENS, AND XIN ZHANG 


7. Numerical illustration 


A basic calibration strategy when presented with asymptotic results such as those given 
in this paper is to assume one can place oneself in the corresponding regime, and then de¬ 
termine model parameters by reading asymptotic coefficient off of market option prices. 
We now illustrate how this strategy can produce positive results, and discuss its limita¬ 
tions, in the context of our Corollary 9 part (i), i.e. when the top of the KL spectrum is 
simple (fi \ = 1). In other words, the idea is to ignore the big O term in the asymptotic 
(52), and calibrate parameters to the remaining coefficients. Denoting the discounted 
log-moneyness log (Soe rT /K) by k for compactness of notation, we thus have, for \k\ suf¬ 
ficiently large. 



(71) 


for two constants L and M, which can, in principle, be read off of market data. By the 
explicit expressions for these constants in (52) in terms of \\ and 8 \, we then express the 
latter in terms of L and M as 


64 T 2 L 4 


(4 — T 2 L 4 ) 2 ' 


4\/2TM\/4 + T 2 L 4 
4 - T 2 L 4 


(72) 


Here we use (71). One notices that, conveniently, \\ can be calibrated using only the 
coefficient L, while given L, 8 \ is then proportional to M. 

At this stage, one may simply conclude that the extreme strike asymptotics given in 
the market are consistent with any Gaussian volatility model whose top of eigenstructure 
is represented by the values computed in the above expressions for Ai and 8 \. However, 
practitioners will prefer to determine a more specific model, perhaps by choosing a classi¬ 
cal parametric one, and using other non-asymptotic-calibration techniques for estimating 
some of its parameters. The expressions in (72) can then be used to pin down other pa¬ 
rameters by calibration, as long as one can relate the model's parameters to its KL pair 
(Ai,^i), whether analytically or numerically. 

We provide illustrations of this strategy in two cases: the stationary Stein-Stein model, 
where the KL expansion is known semi-explicitly, and the Stein-Stein model's long-memory 
version, where the volatility is also known as the fractional Ornstein-Uhlenbeck (fOU) 
model, and the KL expansion is computed numerically. The data we use is also gen¬ 
erated numerically: for each model, we compute option prices and their corresponding 
implied volatilities, by classical Monte-Carlo, given that the underlying pair of stochas¬ 
tic processes is readily simulated. Specifically, in the Stein-Stein (standard OU) case, 10 6 
paths are generated via Euler's method based on discretizing the stochastic differential 
equation satisfied by X started from a r.v. sampled from X's stationary distribution, and 
the explicit expression for log S given X, also approximated via Euler with the same time 
steps; 10 3 time steps are used in [0, T], Option prices are derived by computing average 
payoffs over the 10 6 paths. The details are well known, and are omitted. In the fOU case, 
the exact same methodology is used, except that one must specify the technique used to 
simulate increments of the fBm process which drives X: we used the circulant method. 
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which is based on fBm's spectral properties, and was proposed by A.T. Dieker in a 2002 
thesis: see [1 , 16]. 

Given this simulated data, before embarking on the task of calibrating parameters, to 
ensure that our methodology is relevant in practice, it is important to discuss liquidity is¬ 
sues. It is known that the out-of-the-money call options market is poorly liquid, implying 
that the large strike asymptotics for call and IV prices are typically not visible in the data. 

We concentrate instead on small strike asymptotics. There, depending on the market seg¬ 
ment, options with maturities on the scale of one or several months can be liquid with 
small bid-ask spread for log moneyness as far down as —1. Convincing visual evidence 
of this can be found in Figures 3 and 4 in [20] which report 2011 data for SPX options. 
Based on this, we will use the value —1 as an illustrative lower bound for k: beyond this 
lower bound, liquidity is insufficient to measure IV. This immediately creates a problem 
due to the similar magitude, in this range of k, of the constant term M and the next term in 
the expansion whose order we can only guarantee to be no larger than 1 / \ — k, according 
to the big O term in Corollary 9 part (i). 

In spite of this, calibration of A \ to L as in (72) works well in practice, as our examples 
below show. We also find that the inability to neglect 0^1/ \ Wcj compared to M intro¬ 
duces a bias in the calibration, whose magnitude, while not constant, gives a good idea 
of the relative size of the next term in the small strike expansion, motivating the open 
problem of looking for this next term, and perhaps the following one. 

We begin with the stationary uncorrelated Stein-Stein model with constant mean-reversion 
level m, rate of mean reversion q, and so called vol-vol parameter a. The systems of 
equations needed to perform calibration here have a somewhat triangular structure. Ac¬ 
cording to Section 6, if one is to calibrate q, access to is needed, if one may rely on 
independent knowledge of the level of mean reversion m. Specifically, one would solve 
the following system 


q sin(wT) + wcos(wT) = 0 

Ci fsin(wT) + —(1 — cosiwT))'] = — 

V w / m 

where unfortunately the constant C\ also depends on the unknowns (q, zv) in the follow¬ 
ing non-trivial way: 


1 

S 


|(1 -cos(2a >Tj) + y 


2 w J 2 


sin(2zcT) \ 

2 w ) ' 


Since q is not fixed, the task of determining which value of zv represents the minimal 
solution of the first equation above, given the large number of solutions to the above 
system, is difficult. We did not pursue this avenue further for this reason, and because the 
aforementioned bias due to liquidity means that one would not be in the right asymptotic 
regime. 

Nevertheless, in the Stein-Stein model, our asymptotics allow us to calibrate the vol-vol 
parameter u. The equation for finding a given a measurement of L and prior knowledge 
of q is much simpler. Indeed, since q is assumed given, the base frequency zv is computed 
easily as the smallest positive solution of (56). Independently of this, based on L, the first 
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relation in (72) determine Aj. Then according to equation (61), we obtain immediately 

o 2 = Ai (w\ + q 2 \ . 


We simulated IV data for the call option with m = 0.2 (signifying a typical mean level of 
volatility of 20%), q = 7 (fast mean reversion, every eight weeks or so), and cr — 1.2 (high 
level of volatility uncertainty). How to estimate L from the data is not unambiguous. We 
adopt a least-squares method, on an interval of L-values of fixed length; after experimen¬ 
tation, as a rule of thumb, an interval of length 0.15 provides a good balance between 
providing a local estimate and drawing on enough datapoints. We start the interval as far 
to the left as possible while avoiding any range where the data representing the function 
IhI (, k ) may exhibit clear signs of convexity. This is in an effort to coincide with values 
of k for which sufficient liquidity exists in practice. As a guide to assess this liquidity, we 
use the study reported in [20, Section 5.4], which we mentioned in the introduction. This 
liquidity depends heavily on the option maturity. 

For a one-month (T = 1/12) option, the above calibration method reports a = 1.183, 
based on the interval k G [-0.85,-0.70]. For a three-month option (T = 1/4), the method 
reports cr = 1.193, based on the interval [—1.25, —1.10]. Other calibrations, not reported 
here because of their similarity with these two, show that calibration accuracy increases 
with more liquid options, which is consistent with the intuition that being able to use 
intervals further to the left should allow a better match with the asymptotic regime (71). 
The graphs of the data versus the asymptotic curve, showing excellent agreement for the 


first term La/— /c, are given in Fig. 1 and Fig. 2. 

As mentioned above, this agreement explains why this calibration for cr works well 
even though there is a clear bias between the curve k ha Ly/ —k + M (solid lines in Fig. 
1 and Fig. 2) and the true (simulated) IV data (blue data points in Fig. 1 and Fig. 2). 
One may immediately conjecture that the gap between the solid and discrete lines is rep¬ 
resentative of the term O (1 / y/ —k\ in Corollary 9, in the sense that the next term in the 


expansion might be of the form N / y' —k for some constant N, and that our simulation 
might give us access to the value of N. This does not appear to be quite the case. When 

graphing the expression for I (k) — (k y —k + M j \/ —k using data from the previous 

figure, instead of finding a quantity which is largely constant, we find values which vary 
non-linearly between 0.02 and 0.026, as seen in Fig. 3. This is indicative that, to have 
access to calibration based on the value of M, it might be necessary to find more terms in 
the asymptotic expansion. 

We now explain how to include one more term in formula (3), specifically, how to obtain 
the following expansion: 


+ °(V) (73) 

as k —> — oo when n\ = 1. It is possible to establish formula (73) by using formula (41) 
instead of a weaker two-sided estimate given in (43), and also formula (39) and Theorem 
5.6 in [21], characterizing the asymptotic behavior of log-normal integrals. Combining 
these results, we obtain a sharp asymptotic formula for the density Df with a relative error 
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estimate O ^(logx) 2 J as x —>■ 00 . After that. Theorem 14 in [ '] can be used to establish 

(73). Obtaining the expression for the constant N is tractable but sufficiently non-trivial 
that we omit the details. Even in the case where the volatility process is centered and 
Yi\ = 1, this N depends on all A*. with k > 1 and with k > 2. We also conjecture that 
when n.\ = 1 and 112 = 1, 


I{k) = L\J—k + M + 


N 


+ A+o 


v 7 -fc -k 


.HO 


3/2 


(74) 


as k —> — 00 . This conjecture would not be sufficient to guarantee that including the four 
terms identified above would be enough to provide a good approximation when |fc| is 
restricted to values around 1 or lower, since there a series in negative powers of \J—k will 
diverge unless its coefficients decay fast enough. Yet Fig. 1 and Fig. 3 show, at least in 
those particular examples, that the discrepancy which might correspond to N/V~k + 
P/ (—k) is not large, and that the term P/ (—k) would only be needed to account for 
a moderate uncertainty level on the value of N, ranging from 0.02 to 0.026. This gives 
hope that the first three or four terms in the expansion (74) could provide an excellent 
approximation to the IV even for values of |fc| lower than 1. 

While the conjecture in (74) may be correct, the last analysis above attemps to push 
the boundaries of its possible validity beyond its natural range because of option mar¬ 
ket liquidity considerations. This empirical consideration seems to rely on the class of 
Gaussian volatility models one chooses. This becomes evident in the second class of ex¬ 
amples we present, that of fOU volatility. There, as we will see, the bias between I (k) 
and Ly —k + M cannot reasonably be considered small, so that the conjecture in (74) may 
not be helpful. Nevertheless, our calibration strategy still produces good results when 
applied to the fOU's memory parameter. We now provide the details. 

The continuous-time fOU volatility model was introduced in [8] as a way to model 
long-range dependence in non-linear functionals of stock returns, while preserving the 
uncorrelated semi-martingale structure at the level of returns themselves. This is the 
model for X in (2) where the process Z is a fractional Brownian motion, i.e. the continuous 


= \t — s 


2 H 


Gaussian process started at 0 with covariance determined by E (Zf — Z s 

with "Hurst" parameter H G (0.5,1). In [7], it was shown empirically that standard statis¬ 
tical methods for long-memory data are inadequate for estimating H. This difficulty can 
be attributed to the fact that the volatility process X can have non-stationary increments. 
In addition, some of the classical methods use path regularity or self-similarity as a proxy 
for long memory, which cannot be exploited in practice since there is a lower limit to 
how frequenty observations can be made without running into microstructure noise. To 
make matters worse, the process X is not directly observed; in such a partial observation 
case, a general theoretical result was given in [32], by which the rate of convergence of 
any estimator of H cannot exceed an optimal H-dependent rate which is always slower 
than N _1//4 , where N is the number of observations. Given the non-stationarity of the pa¬ 
rameter H on a monthly scale, a realistic time series at the highest observation frequency 
where microstructure noise can be ingored (e.g. one stock observation every 5 minutes) 
would not permit even the optimal estimators described in [32] from pinning down a 
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value of H with any acceptable confidence level. The work in [7] proposes a calibration 
technique based on a straightforward comparison of simulated and market option prices 
to determine H. Our strategy herein is similar, but based on implied volatility 
We consider the fOU model described above with the following parameters: T = 
1/12, m = 0.2, q = 7, cr = 1.2, with different values of the Hurst parameter H £ 
{0.51,0.55,0.60,0.70,0.80}. As mentioned above, we simulated option prices using stan¬ 
dard Monte Carlo, where the fOU process is produced by A.T. Dieker's circulant method. 
Since the values of Aj for each H > 0.5 are not known explicitly or semi-explicitly, we 
resorted to the method developed in by S. Corlay in [9] for optimal quantification: there, 
the infinite-dimensional eigenvalue problem is converted to a matrix eigenvalue problem 
which uses a low-order quadrature rule for approximating integrals (a trapezoidal rule 
is recommended), after which a Richardson-Romberg extrapolation is used to improve 
accuracy. We repeat this procedure for the fOU process with the above parameters, for 
each value of H from 0.50 to 0.99, with increments of 0.01. The corresponding values we 
obtain for Ai in each case are collected in the following table: 


H = 

0.50 

0.51 

0.52 

0.53 

0.54 

0.55 

0.56 

0.57 

0.58 

0.59 

Ai = 

0.00713 

0.00697 

0.00682 

0.00667 

0.00653 

0.00638 

0.00624 

0.00610 

0.00597 

0.00584 

H = 

0.60 

0.61 

0.62 

0.63 

0.64 

0.65 

0.66 

0.67 

0.68 

0.69 

Ai = 

0.00571 

0.00558 

0.00546 

0.00534 

0.00522 

0.00511 

0.00499 

0.00488 

0.00478 

0.00467 

H = 

0.70 

0.71 

0.72 

0.73 

0.74 

0.75 

0.76 

0.77 

0.78 

0.79 

Ai = 

0.00457 

0.00447 

0.00437 

0.00428 

0.00419 

0.00410 

0.00401 

0.00392 

0.00384 

0.00376 

H = 

0.80 

0.81 

0.82 

0.83 

0.84 

0.85 

0.86 

0.87 

0.88 

0.89 

Ai = 

0.00368 

0.00360 

0.00353 

0.00345 

0.00338 

0.00331 

0.00324 

0.00317 

0.00311 

0.00305 

H = 

0.90 

0.91 

0.92 

0.93 

0.94 

0.95 

0.96 

0.97 

0.98 

0.99 

Ai = 

0.00299 

0.00293 

0.00287 

0.00281 

0.00275 

0.00270 

0.00265 

0.00260 

0.00255 

0.00250 


Our illustration of the calibration method then consists of starting with simulated IV 
data for a fOU model with fixed H, then implementing the same procedure as in the Stein- 
Stein example to determine the value of \\ by using an interval of length 0.15 which is 
realistic in terms of liquidity constraints, and then matching that value of \\ to the closest 
value in the above table, thereby concluding that the simulated data is consistent with 
the corresponding value of H in the table. The results of this method are summarized as 
follows 


True H 

0.51 

0.55 

0.60 

0.70 

0.80 

Interval used 

-1.0,-0.85] 

-1.0,-0.85] 

-1.0,-0.85] 

-1.0,-0.85] 

-1.0,-0.85] 

measured Aj 

Ai = 0.0071 

Ai = 0.0064 

Ai = 0.0055 

Ai = 0.0045 

A x = 0.0037 

calibrated H 

0.50 

0.55 

0.62 

0.71 

0.80 


Our method shows a good level of accuracy, though it does not appear capable of dis¬ 
tinguishing between the case of no memory H = 0.50 and a simulation with very low 
memory (H = 0.51). The accuracy of the method has some sensitivity to the intervals 
being used, particularly for small true values of H: in that case the first term in the im¬ 
plied volatility asymptotic curve seems to dictate a need for including data points for log- 
moneyness k near —1. For higher values of H, we noted greater robustness to the interval 
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being used to estimate L. For example, with a true H = 0.70 used in simulation, cali¬ 
brated values of H = 0.69 and H = 0.70 were obtained with the intervals [—0.87, —0.72] 
and -0.90,-0.75] respectively. 

More generally, one notes that the bias between the curve La/ —k + M and the simulated 
IV data is much larger when H > 0.5 than in the classical Stein-Stein case where H = 0.5, 
as illustrated in Fig. 4. This provides an even stronger call than in the case H = 0.5 to 
compute the constants N and P in the conjecture in relation (74), as functions of the first 
few KL eigen-elements. 

We also computed the fOU's KL eigenvalues A 2 and A 3 for each H E {0.50,0.51,..., 0.99}, 
though these are not used in our calibration. They are nevertheless instructive since they 
indicate at what speed the KL expansions of X and of the integrated variance Y might con¬ 
verge. Without reporting all the values, we find that the ratio A 2 /A 1 decreases from 0.11 
for H = 0.50 to 0.0010 for H = 0.99, while for A 3 /A 1 , the values are even smaller, decreas¬ 
ing from 0.033 to 0.00012 over the same range of H. If the values of N/L and P/L in our 
conjecture (74) are of the same order as the ratio A 2 / Ai, this could give additional hope 
that computing N, P in the fOU case might be enough to implement a two-parameter 
calibration, with only a few values of A, needed to approximate M, N, P. 

Finally, by comparing Fig. 1 with the first graph in Fig. 4, where all OU parameters 
remain the same except for H which goes from 0.5 to 0.51, one notes that the difference 
I (. k ) — (La /—k + M) jumps from values on the order of 0.02 to much larger values, on 
the order of 0.2. This indicates that the constant N in formula (73) might have a discon¬ 
tinuous dependence on H as H approches 0.5; this phase transition between the case of 
no memory (H = 0.5) to that of long memory (H > 0.5) is presumably due to a jump 
in the value of KL elements as H approaches 0.5. This could lead to a method for deter¬ 
mining whether long memory is present, but in the case of IV asymptotics it would be 
constrained by liquidity considerations. We will investigate all these questions in another 
article. 
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FIGURE 1. One-month IV for Stein-Stein model with parameters m = 0.2, 
q = 7,a = 1.2 



FIGURE 2. Three-month IV for Stein-Stein model with parameters m = 0.2, 
q = 7,a = 1.2 



FIGURE 3. Graph of I (k) — (L\/—k +Mj \/—k with simulated data 
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Figure 4. IV for fOU model with H = 0.51, H = 0.55, H = 0.60, H = 0.70, 

H = 0.80 
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